Appeasing the Phantom Menace? 
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An induced gravity brane- world model is considered herein. A Gauss-Bonnet term is provided 
for the bulk, whereas phantom matter is present on the brane. It is shown that a combination of 
infra-red and ultra-violet modifications to general relativity replaces a big rip singularity: A sudden 
singularity emerges instead. Using current observational data, we also determine a range of values 
for the cosmic time corresponding to the sudden singularity occurrence. 
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I. INTRODUCTION 

From current observational data [3-[3l , it is now widely 
accepted that the universe is undergoing a state of accel- 
eration. The simplest setup to describe this acceleration 
is by means of a cosmological constant, with an equa- 
tion of state pa = —p\, where p\ is the pressure and 
Pa is the energy density of such a cosmological constant. 
There are, however, other candidates that astronomical 
observations still allow. Namely, a dark energy compo- 
nent 0, [H whose equation of state parameter, w, is very 
close to —1, where w is the ratio between the pressure of 
dark energy, p, and its energy density p. The relevance of 
this point is that dark energy can lead to quite different 
scenarios concerning the future of the universe. To be 
more precise: if w > —1, dark energy corresponds to a 
quintessence fluid Q;iftu = — 1, itisa cosmological con- 
stant and the universe would be asymptotically de Sitter; 
if w < —1, dark energy corresponds to a "phantom" con- 
tent d). 

In the context of dark energy cosmology, the study 
of gravitational theories that entail spacetime singular- 
ities at late-times has made a considerable progress in 
the last years. Being more specific, the following cases 
have been established within a Friedmann-Lemaitre- 
Robertson- Walker (FLRW) framework 1 Q : 

• Big rip singularity - a singularity at a finite cosmic 
time where the scale factor, the Hubble rate and 
its cosmic time derivative diverge, (til [9l Hoj|; 

• Sudden singularity - a singularity at a finite scale 
factor and in a finite cosmic time where the Hubble 
rate is finite but its cosmic time derivative diverges 
EL El; 
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1 For an alternative classification of the dark energy related singu- 
larities see Ref. [jj . 



• Big freeze singularity - a singularity at a finite scale 
factor and in a finite cosmic time where the Hubble 
rate and its cosmic time derivative diverge fl3l j: 

• Type IV singularity 2 - a singularity at a finite scale 
factor and in a finite cosmic time where the Hubble 
rate and its cosmic derivative are finite but higher 
derivative of the Hubble rate diverges. These sin- 
gularities can appear in the framework of modified 
theories of gravity . 

Subsequently, it has become of interest to determine 
under which conditions any of the above singularities 
can be removed, or at least appeased in some manner. 
A significant range of approaches and corresponding re- 
sults have been contributed to the literature |14l - l2ll | . Our 
purpose in this paper is to participate in that line of 
investigation, namely addressing the following question: 
Since the big rip singularity occurs at high energy and 
in the future, could we expect that a combination of IR 
and UV effects removes or replaces it? To this aim, we 
employ in this paper a Dvali-Gabadadze-Porrati brane- 
world (DGP) model [H,[23[ (see also [3), that includes 
a phantom matter fluid (that emulates the dark energy 
dynamics), where the five-dimensional bulk is character- 
ized by a Gauss-Bonnet (GB) term. Generally speak- 
ing, within a DGP brane component, we have a setting 
where infra-red (IR) effects modify general relativity at 
late-time, whereas with the GB ingredient, ultra-violet 
effects are present for high-energy scales |25|-|27| . 

Before proceeding into a more technical discussion, it 
could be of interest to further add the following about 
having a phantom matter fluid for the brane. On the 
one hand, dark energy component with w < — 1, i.e. a 
phantom energy component has not yet been excluded 
by the recent results of WMAP5. For example, the 
WMAP5 data (in combination with other data) for a 
standard FLRW universe with spatially flat sections, 
filled with cold dark matter and a dark energy component 
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with a constant equation of state parameter, w, predicts 
-0.097 < 1 + w < 0.142 for CMB, SNIa and BAO data 
which gives the most stringent limit, while WMAP data 
alone predicts for this model —0.78 < 1 + w < 0.68. 
For more details, see [28|. On the other hand, to inves- 
tigate future singularities, a perfect fluid is satisfactory 
and therefore we have not given an explicit action for the 
phantom matter in terms of a minimally coupled scalar 
field (with the opposite kinetical term) or through more 
general scalar field actions like a k-essence action (see 
Section II) . Let us also add that it is well known that the 
DGP brane has a branch with an unstable mode solution 
(i.e., a 'ghost') |29l |. It may therefore be questionable why 
to initiate a study within a DGP setting or even insert 
phantom matter 3 . Our point is that, in spite of these 
open lines, the features characterizing the DGP as well 
as GB elements can provide a framework to investigate 
the intertwining of late time dynamics and high energy 
effects. An interesting ground to test it is with a phan- 
tom fluid, since this matter in a standard FLRW setting 
induces a big rip singularity. Eventually, the unresolved 
issues for the DGP brane will be eliminated and results 
such as the one we bring here will increase in interest. 

This paper is then organized as follows. In Sect. II, 
we describe our DGP-GB model filled with matter and 
a phantom fluid and present the analytical solutions. In 
Sect. Ill, we analyze how the big rip is replaced by a sud- 
den singularity. Finally, we summarize our work in Sect. 
IV, presenting also some possible lines to subsequently 
investigate. 



will define their sum as p m \ i.e. p m = pb + Pcdm- On the 
other hand, we will consider dark energy to correspond 
to phantom energy. Finally, the total energy density on 
the brane can be written as 4 



PmO 



Pda 



a 3(w+l) > 



(3) 



where w, p m o-, Pdo are constants and 1 + w < 0. 

It should be noted that from ([1]), we can obtain the 
known self- accelerating DGP solution [23j (+ sign in 
Eq. ((T|) with a = 0), while the normal branch is retrieved 
for the — sign with a = 0; cf. [27j for more details and 
notation. 

Let us then address equation (p} analytically, selecting 
the + sign, adopting the approach presented in [27[ • This 
equation with the matter content ([3]) can be expressed as 



e 2 (z) = n m (i + z) 3 + n d (i + z) 3( - 1+ ^ 
+2^n7 c [i + n a E 2 {z)} E(z), 

where E(z) = H/H n , z is the redshift and 
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Q a = -aH a . 



(5) 



Evaluating the Friedmann equation @ at z — gives 
a constraint on the cosmological parameter of the model 



i = a r 



(0) 



II. THE DGP-GB MODEL WITH PHANTOM 
MATTER 

The generalized Friedmann equation of a spatially flat 
DGP brane with a GB term in a Minkowski bulk can be 
written as [25[ (see also [26^ 



H2 = fp ± ~ 

Qr c r c 



-aH 2 



H. 



(1) 



where a is the GB parameter, k 2 is the five-dimensional 
gravitational constant and r c is the crossover scale. 
Moreover, p stands for the total energy density of the 
brane. Therefore, for a late-time evolving brane the en- 
ergy density is well described by 



P = Pb+ Pcdm + , 



(2) 



where pb, p c dm and pd corresponds to the energy density 
of baryons, cold dark matter and dark energy, respec- 
tively. As pb and p a dm are both proportional to a -3 , we 



For f2 Q = 0, we recover the constraint in the DGP model 
without UV corrections. Coming back to our model, if 
we assume the dimcnsionless crossover factor Q rc to be 
the same as in the self-acceleration DGP model, then the 
similarities with a spatially open universe are made more 
significant from the GB effect, since fl a > 0. 

In order to obtain the evolution of the Hubble rate as 
a function of the total energy density of the brane, we 
introduce the following dimcnsionless variables: 



H 



8 a 
3 r c 



H = 2n a ^/aT c E{z), 



(7) 



32 n 2 a 2 
27~^r P 



n d (i + z) 3 ^ + n m (i + z) 3 } , (8) 



3r 



(9) 



In terms of these variables, the modified Friedmann equa- 
tion then reads: 



H 3 - H 2 + bH + p = 0. 



(10) 



For a theoretical criticism on phantom energy models, see how- 
ever, 1301. 



4 Quantities with the subscript denote their values as observed 
today. 
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Notice that there is a change of sign with respect to equa- 
tion (17) in (27}- The number of real roots is determined 
by the sign of the discriminant function ./V defined as 5 



N = Q 3 +R 2 , 



(11) 



where Q and R are, 



R = -l b -y + Tr (12) 



For the analysis of the number of physical solutions of 
the modified Friedmann equation (|10[) . it is helpful to 
rewrite N as 



N = ^(p-piXp-pa), 



(13) 



where 



Pi = -^{b-f[l + V(l-36) 3 ]}, (14) 



a unique solution for the cubic Friedmann equation 
(|10p , because N is a positive function for this case. 
So the solution reads, 



Hi = - 
where r\ is defined by 



1 



cosh(?7) 



■R 



sinh(r/) 



N 



(16) 



(17) 



and 77 > 0. When 77 — > 0, the energy density of the 
brane approaches pi. This solution has a negative 
Hubble rate and therefore it is unphysical for late- 
time cosmology. 

• During the limiting regime, p — pi, the function N 
vanishes, and there are two real solutions 



H x = -- (2VI - 36- l) 



(18) 



P2 = -i{fe-^[l-v / (l-36) 3 ]}. 



(15) 



Hence, if TV is positive then there is a unique real solution. 
If N is negative, there are three real solutions, and finally, 
if N vanishes, all roots are real and at least two are equal. 

An approximated bound for the value of 6 can be es- 
tablished noticing that 6 is proportional to Cl rc through 
(|9]). Hence, from the equivalent quantity for the £l rc in 
the DGP scenario for the self-accelerating branch [32[ and 
the constraint on the curvature of the universe 6 (see , 
e.g.), its value should be small. These physical solutions 
can be included on the set of mathematical solutions with 
< b < j. Therefore, for the remaining of this letter, 
we shall study in detail this setup. For completeness the 
other cases are summarized in Table HI 

For < b < j the values of pi and p2 in equations (fT4|) 
and (flip)) are real. More precisely, in this case pi > and 
p2 < 0. The number of solutions of the cubic Friedmann 
equation (|10[) will depend on the values of the energy 
density with respect to pi. As the (standard) energy 
density redshifts backward in time (i.e. it grows), we 
can distinguish three regimes: (i) high energy regime: 
pi < p, (ii) limiting regime: p = pi, (iii) low energy 
regime: p < pi. 

• During the high energy regime, the energy density 
of the brane is bounded from below by pi . There is 



The modified Friedmann equation J 1 1 ) will be solved analyti- 
cally following the method introduced in Ref. [2^1 for the normal 
DGP-GB branch. For a semi-analytical approach see the second 
reference in Ref. [25l |. 

At this respect, notice that at z = 0, the term 2^/f2r e (1 + O a ) 
mimics a curvature term on the modified Friedmann equation. 



H 2 = -(VT=3b+l). 



(19) 



The solution Hi is negative and so it is also not 
relevant physically. 

• For the low energy regime, p < pi. In this case 
the function N is negative, and there are three dif- 
ferent solutions. One of these solutions is negative 
and corresponds to a contracting brane, while the 
other two positive solutions correspond to expand- 
ing brancs. Let us be more concrete: 

The solution that describes the contracting brane 
is similar to the corresponding solution of the high 
energy regime: 



Hi = -i 
3 



where 



cos(0) 



-R 

7^ 



2VT - 36 cos [ - I - 1 



sin(0) 




(20) 



(21) 



and < 9 < 9q. The parameter 9 = is defined as 
in equation (|2ip in which the value of p reaches pi , 
and the parameter 9q corresponds to p = 0. For this 
solution Hi is negative and hence not suitable for 
the late-time cosmology. In addition, the solution 
approaches the same Hubble rate at 9 = as the 
limiting solution ([T8]) . 

On the other hand, the two expanding branches are 
described by 



2\/l - 36 cos 



1 



(22) 
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b 


pi and p2 


Solutions for H 


rj and i? 


Description 


4 < b < 3 


Pi < 0, P2 < 


iii = — vl — 3ocosn(j) — 1] 


i / \ R 

cosh(a) = , 

sinh(a) = « ^ks, 
qo < a, where 


Hi < 0; contracting 
brane. 




pi = p2 = 


fli = -i[(l + 27p)-*-l] 




Hi < 0; contracting 
brane. 


\<b 


pi , p2 : complex 
conjugates 


Hi = -i[2V36- lsinh(f) - 1] 


cosh(tf) EE ^ 
sinh(tf) = 

< $, where 


Hi < 0; contracting 
brane. 



TABLE I: Solutions for the algebraic equation with different ranges for b; See also Eqs. (|16[) - (|23|l . 




FIG. f : Plot of the redshift at which the total energy density of the brane reaches its minimum value. 



H* = - 



cos 



(23) 



For 0—^0 the energy density p approaches p\ 
and the low energy regime connects at the limiting 
regime with the solution (fl~9| where both solutions 
coincide; It can be further shown that Hi < H3. 
The more interesting solution for us is the brane 
expanding solution described by Eq. (f2"2")l : This 
solution constitutes a generalization of the self- 
accelerating DGP solution with GB effects but 
within our herein 7 model. 



III. THE BRANE ESCAPING THE BIG RIP 
AND PLUNGING INTO A SUDDEN 
SINGULARITY 

The usual self-accelerating DGP-GB solution is known 
to have a sudden singularity in the past, when the brane 
is filled by standard matter [25| . Herein, we consider in- 
stead the brane filled with CDM plus phantom energy, 



performing an analysis aiming at future cosmic time. 
Starting from the modified Friedmann equation (|10[) , it 
can be shown that the first derivative of the Hubble rate 
is 



H 



3[2H- £(1 + 8aH 2 )}' 



(24) 



where p is given by the equation of the energy conserva- 
tion 



p + 3H[p m + (1 + w)p d ] = 0, 



(25) 



for the total energy density of the brane. Therefore, p 
decreases initially until the redshift reaches the value z*, 



1 



-(1 + tu) 



fid 

- - 1 1 



(26) 



where p — 0. Afterwards, the phantom matter starts 
dominating the expansion of the brane; indeed the brane 
starts super-accelerating 8 (0 < H) and the total energy 



7 It should be clear that we have herein on the paper a DGP-GB 
setting with phantom matter and therefore our brane does not 
has, technically speaking, a self-accelerating phase, asymptoti- 
cally approaching a de Sitter stage in the future. 



8 Notice that the denominator of equation ilL'tt is always positive 
for the self-accelerating DGP-GB solution This can be 

proven by using Eqs. 1)240 . 122t and the range of the variable 9. 
Indeed, it can be shown that < 3[2H - i (1 + SaH 2 )] < i=^. 

This means that in our model H is negative for z* < z and 
positive for z < 2*. 
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density of the brane starts growing as the brane expands. 
In figure [TJ we plot the redshift z+ versus the observa- 
tional values of the parameters 9 fid, flm and w. As can 
be noticed in this figure, the total energy density of the 
brane would start increasing only in the future. Further- 
more, the larger is \w\, the sooner the phantom matter 
would dominate the expansion of the brane. Q m has the 
opposite effect while fl d has a much milder effect on z+. 



the first derivative 



Substituting Eq. fl25j) into Eq. (|2l 
of the Hubble rate reads, 
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k\H[{1 + w)p d + p m 
2H - ±{l + 8aH 2 ) 



(27) 



This equation shows that when the Hubble rate ap- 
proaches the constant value 10 , 



8a 



(28) 



the first derivative of the Hubble parameter, H, diverges, 
while the energy density of the brane remains finite. 
Thus, instead of a scenario where the energy density on 
the brane blucshifts and eventually diverges, with a big 
rip singularity emerging, we find that (i) a finite value 
of the (dimensionless) energy density p — pi (the same 
applies for the pressure) , and (ii) a finite value for the (di- 
mensionless) Hubble parameter H = H 2. Those values 
are reached in the limiting regime described by Eq. (|19p . 
And (iii) the first derivative of the Hubble parameter di- 
verges at that point. Therefore, the energy density is 
bounded, i.e. the limit z — > — I or a — > 00 cannot be 
reached, and instead of a big rip singularity we get a 
sudden singularity, despite the brane being filled with 
phantom matter. 



9 We estimate in figure [Tl by choosing values for Sl m , f2,j and 
10 in accordance with the WMAP data . We know that those 
parameters are model dependent and therefore a best fit analysis 
of the brane expansion using the currently available observational 
data (for example SNIa, BAO and CMB) would provide a much 
better analysis of z+. However, this analysis is far beyond the 
scope of the current work. 
10 The Hubble rate given in Eq. 1128ft can be mapped into the di- 
mensionless Hubble parameter given in Eq. Ill9t for the limit- 
ing regime, reached at the constant dimensionless energy density 
P = P1- 



FIG. 2: Values for the redshift at the sudden singularity, z S mg, 
for fixed value of Q d = 0.75, fi m = 0.24 and fl rc = 5.95 x 10~ 6 
with respect to the parameter w. 

The results above can be further elaborated, by estab- 
lishing when the sudden singularity will happen (i.e., at 
which redshift and cosmic time values): 

• In order to obtain the redshift z s ; ng where the sud- 
den singularity takes place, we equate the dimen- 
sionless total energy density of the brane (JSJ to its 
values at the sudden singularity (cf. Eq. (fl4| ). The 
allowed values for z sulg are plotted in Fig. [2] for the 
various values of the equation of state parameter 
w. 

• Finally, using the relation between the scale factor 
and the redshift parameter, a(t) = 1/(1 + 2), one 
can write the Hubble rate as a function of the red- 
shift parameter and its cosmic time derivative as 
follows 



H = - = , 

a 1 + z 



(29) 



and then integrating this equation, the cosmic time 
remaining before the brane hits the sudden singu- 
larity reads 



(ising — to)Ho 



dz 



{l+z)E{z) 



(30) 



In the previous equation to and i s i ng indicate the 
present time and the time at the sudden singular- 
ity, respectively. We can plot i sing for fixed values 
of (f2 rn , fid, ^a); See Fig. [3] As we can notice from 
Fig. 02 the closer is the equation of state of the 
phantom matter to that of a cosmological constant, 
the farther would be the sudden singularity. The 
same plot is quite enlightening as we can compare 
the age of the universe, essentially Hq 1 , to the time 
left for such a sudden singularity to take place on 
the future of the brane. As we can see, such a sud- 
den singularity would take place roughly in about 
0.1 Gyr. 
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It is of interest to compare this estimate with the time 
at which a big rip would take place in a standard four- 
dimensional universe filled with phantom matter whose 
equation of state is constant. In ref. [Toj . the time re- 
maining for such a big rip singularity was estimated to be 
about 22 Gyr. However, notice that in [Toj . w = —1.5 and 
Q m = 0.3 was used, whereas we have w € [—1.1, — 1-01], 
Q m = 0.24; Our setting is also different. Within the 
framework of [l(|, for w = —1.1 and O m = 0.24, one 
retrieves a time of about 80 Gyr. In comparison, for 
the model presented in this paper, with w = —1.5 and 
O m = 0.3, Q<j = 0.68, il a = 0.005, the singularity 
emerges in 0.2 Gyr. On the whole, although the universe 
would meet a "less severe" singularity in the future, this 
would happen much sooner that for the big rip case. In 
fioj . there are a few predictions concerning possible as- 
tronomical events that would indicate the emergence of 
the singularity herein exposed. Within our model, some 
of those events or other of similar impact could occur 
rather earlier. 




-1.1 -1.08 -1.06 -1.04 -1.02 

W 



FIG. 3: Variation of the cosmic time left before the brane 
hits the sudden singularity, (f s i ng — to)Ho (in dimensionless 
units), for fixed values of fid = 0.75, Q m = 0.24 and Q, Tc = 
5.95 x 10 -6 with respect to the equation of state parameter 
w. 



IV. CONCLUSIONS AND OUTLOOK 

The study of singularities occurring in the future is a 
subject of interest and several proposals for their removal 



or substitution have been advanced [lj-|2l|. In this pa- 
per, we have considered a specific case and investigated 
whether a composition of specific IR and UV effects could 
alter a big rip singularity setting. More concretely, we 
employed a simple model: a DGP brane model, with 
phantom matter and a GB term for the bulk. The DGP 
brane configuration has relevant IR effects, whereas the 
GB component is important for high energies; phantom 
matter in a standard FLRW model is known to induce 
the emergence of big rip singularities [f| . 

Our analysis indicates that the big rip can be replaced 
by a sudden future singularity, through some intertwin- 
ing between late-time dynamics and high energy effects. 
Subsequently, we determine values of the redshift and 
cosmic time, before the brane reaches the sudden singu- 
larity. These results can be contrasted with those, e.g., 
in [l(| for the big rip occurrence in a FLRW setting. 

We are aware that the herein conclusions are based 
in a rather particular result, that was extracted from 
a specific model. Subsequent research work would as- 
sist in clarifying some remaining issues. For example, 
it would be interesting to further study if and how other 
singularities can be appeased or removed by means of the 
herein combined IR and UV effects. On the other hand, 
it might be interesting to consider a modified Hilbcrt- 
Einstein action on the brane, which in addition could al- 
leviate the ghost problem present on the self-accelerating 
DGP model by self-accelerating the normal DGP branch 
HU , and see if some of the dark energy singularities can 
be removed or at least appeased in this setup. We hope 
to report on these lines in a forthcoming publication. 
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